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ABSTRACT 

The effect of various jamming waveforms on conventional 
binary incoherent digital receivers was analyzed in terms of 
resulting receiver performance (i.e., receiver probability 
of error) . 

Probability density functions associated with the test 
statistic generated by incoherent receivers under the influ- 
ence of noise and jamming have been obtained. 

Due to the complexity of the mathematical expressions 
specifying receiver probability of error in closed form, no 
attempt has been made to obtain absolute optimum jamming 
waveforms operating against binary incoherent receivers. 
Therefore near optimum jammer signals were proposed, 
studied, and evaluated. 

The effect of a varying threshold on receiver perform- 
ance was investigated and a jamming strategy involving use 
of an FM jammer was considered, and its effect evaluated. 

Graphical results are presented that highlight the math- 
ematical results obtained. 
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I. INTRODUCTION 



A. COHERENT CORRELATION RECEIVER 

The optimum receiver which will detect the presence of a 
signal or discriminate between two different signals in the 
presence of additive white Gaussian noise, is the well-known 
coherent correlator receiver structure which has two equiva- 
lent forms, as shown in Figure 6.1 [Ref. 2]. 

The binary communication problem is modeled using 
hypothesis testing theoretic principles in which one of two 
signals, S, (t) or S 0 (t), is received in the time interval 
[0,T]. Due to the presence of the noise, the observable 
signal r(t) takes on one of the two following forms 



where the noise n(t) is assumed to be a sample function of a 
stationary white Gaussian process having zero mean and power 
spectral density (PSD) level N 0 /2 watts/Hz. 

The optimum receiver (in the sense of minimum prob- 
ability error, Pe) generates the test statistic G, where G 
is given by 



H, : r(t) = S, (t) + n(t) 



0{t4T 



or 



H, : r(t) = S„(t) + n(t) 



0 < t £T 




or equivalently 



G = <r,S,> - <r,S.) * jfilS.f- ||S,f] 
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and compares it to a threshold in order to decide whether 
signal S i ( t ) or S e (t) has been transmitted in the time 
interval [0,T]. The inner product notation (♦,*) implies 



(x,y) = fx(t)y(t) dt 
Jo 



and the norm notation ||*|| indicates 
: (x,x) = fx*( t) dt 



x 



In order to compute receiver probability of error, it is 
necessary to determine the probability density function of G 
conditioned on the hypotheses H, and H 0 . Denoting these 
density functions P ( (G) and P 0 (G) depending on whether S ( (t) 
or S 0 (t) respectively was transmitted, the error probability 
P e can be expressed as 



n 



P e = P(H , ) 



r°° 



P, (G) dG + P(HJ 



P„ (G) dG 



dec 



J r 



where 



*= ik 

' 2 . 



L 



P(H q ) 
P(H , ) 



and P(H|) denotes the priori probability that the hypothesis 
H | is true. It can readily be shown that if P(H,) = 1/2, 



CD 



where 



^ |V27T 

[0-p)e/nJ 

J_ fL, 



-Vz , 

^ d; 



(i.i) 



2 



j[ sf'Ct) * S„ 2 (t)j dt 



and 



pt 



£ / S, (t)S 0 (t) dt 
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We can interpret E as the average energy per bit and p as 
the normalized signal cross correlation. Using the 

following definition of the error function 

f‘ * | 

erfc * (x) "lw e 

Equation 1.1 can be written as 



? e - erfc* 



Jd -p)E/\] a 



From the above equation it can be noted that as (l-p)E/N 0 
increases, the error probability decreases. For fixed E/N 0 
therefore, the optimum signals choice is one for which P = 
-1. This occurs whenever S,(t) = -S 0 (t). A system employing 
this choice of signal is known as an optimum antipodal 
signaling binary communication system. 



B. JAMMING OF COHERENT RECEIVERS 
1 . Gen e ral 

The objective of this section is to summarize previously 
derived results on the performance of the coherent binary 
(optimum) receiver of Figure 6.1 operating in a jamming 
environment. That is, the transmitted signal is interfered 
by the presence of a jamming waveform as well as additive 
white Gaussian noise [Ref. 1]. The signal appearing at the 
front-end of the receiver can be mathematically modeled by 

H, : r(t) = S,(t) + n(t) + n.(t) 

or 

H 0 : r(t) = S 0 (t) + n(t) + n.(t) 

where n.(t) is a jammer waveform modeled as deterministic, 
\) 

yet unknown to the receiver. 



0 4t ^T 



0 ^ t 4 t 



( 1 . 2 ) 
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2 . Jammer Optim izat ion 

The effect of a jammer waveform on the optimum cohe- 
rent receiver is now presented in terms of its impact on the 
receiver probability of error. 

In order to determine the performance of the 
receiver, we analyze the test statistic G, which is given by 

G = (r,S d ) 4[N| ~ || S, | ] 



with r(t) now given by Equation 1.2 . Due to the assumption 
of a Gaussian zero mean noise having PSD level N c /2 
watts/Hz, it is simple to show that G is a conditionally 
Gaussian random variable with mean value 

E (g|S.| = (S-,Sj) * (n.,Sj) 4[||s/ ' ll S 'l| 2 ] • 1 = 0 ’ 1 
and variance 

Var{ G |Si j ^IMI* 1 = O. 1 

With the mean and variance known, the resulting conditional 
density function of G is given by 




From the general probability 
1.1) and an approriate change 
when P (H , ) = P(H 0 ) = 1/2 the 



r 




of error expression (Equation 
of variables in the integrals, 
expression for Pe becomes 




-co 
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where 



Vz. 



No S 



M - > ISdl|A 

N d - “ 2 ~ ; d = (n.,Sj) 



From the jammer standpoint, the optimum jammer waveform must 
be chosen in such a way that it maximizes the receiver prob- 
ability of error. A determination of the optimum jammer 
waveform can be made by first evaluating the derivative of 
Pe with respect to d, which is the cross correlation between 
the jammer waveform and the signal difference Sj(t). Based 
on Equation 1.3, we can show that 



3Pe_ 




sinh SyNj • d 



3d JT TT 

Due to the behavior of the hyperbolic sine function, we 
observe that 



BRe . 
3d " 



>0 

= 0 

<0 



, d > 0 
, d = 0 

, d <0 



and 



3P e 

3d 1 



JTtT 



exp [ - S* (d + N d f/2 ] > 0 



for all values of d. Since Pe is. a monotonic function of d, 
it is apparent that making d as large as possible in magni- 
tude results in the largest possible increase in Pe . In the 
limit as | d | — > oo , it is easily seen that Pe --> 1/2. 
However, from the Cauchy - Schwarz inequality 
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(1.4) 




with equality if n . ( t ) is proportional to S^(t). 

I In. I I = / R. , where P*. is the jammer energy, 

j * 1 j j 

| d | - - > oo implies that P n . - - > oo when | | Sj | | < co 

■J 

it is not possible to have infinite jammer energy, 

constrain n.(t) such that P^. is finite, 
u o 

Cauchy- Schwarz inequality, Equation 1.4 can be mad 
equality if 



Defining 
the term 
•Since 
one must 
From the 
e into an 



n.(t) = k*S d (t) 



where K is a constant of proportionality. 
K must be set to the value J P^. / | | Sj | | . 
ized by setting 



If | |n. | | = /p7, 

'J j 

Thus |d| is maxim- 



(1.5) 

n j' (t) =/^T IISjII 



and from the above discussion, this results in Pe being 
maximized also. Substituting Equation 1.5 into the prob- 
ability of error expression (Equation 1.3) yields 




It is possible to put Equation 1.6 in a more meaningful form 
by defining 

E * 

N,> “ : signal to noise ratio 
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jammer to signal ratio 



and observing that 
. 2 . 



|| s d|| = I [S, (t) - S 0 (t )] dt = 2 E (1 -p) 

J O 

Then, the probability of error (Equation 1.6) becomes 






■co 



-JSN R(Jl-P t J2 JSr) 



- x i /2 



L 



/2 7C @ dx +j 

JjSNRC/^p -/2JSfr ) - 



/27f <g 



- x*/2 



dx 



(1.7) 



Analysis of Equation 1.7 indicates the fact that for a value 
of JSR beyond a so-called 'break point', that is JSR > 
(l-p)/2, an increasing SNR causes Pe to become worse, i.e., 
it increases. 

3. O ptimum Jammer Wave f orm s for PSK, FSK and ASK 

The effectiveness of a deterministic jammer waveform on 
various modulation techniques will now be presented for PSK 
modulation, where 

S,(t) = A sinW c t ; S 0 (t) = A sin(W c t + 7H) O^t^T 

with the constraint that WcT = n 70 ; n an integer. Also, 
results for FSK modulation will be shown, where 



S,(t) = A sin W,t ; S 0 (t) = A sin W a t 



0 ^T 
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with the constraint that (W| - W 0 )T = 1 7H ; (W t 
m7L ; 1 and m integers. Finally for ASK modulation 

S,(t)=A,S(t) ; S 0 (t) = A c S(t) 0 A t AT 



where we assume that | |S| | < cO and for convenience 

> A c . The optimum jammer waveforms and resulting 
performances can be obtained using Equations 1.5 
Thus the optimum jammer for PSK is given by 



n.(t) = J P_ 



cos W c t 0 At £ T 



j ’ Vt 

Since P = -1 for PSK, we have 



P = — 
r e 2 



■CD 



r 



1 x 2 /2 

J27C^ 



dx + 



J^srn (i -Jj sr) 



4 



2 SNR (1 + 7j S R) 






J27L 



dx 






For FSK modulation, the optimum jammer is given 

2 



n.(t) = 

j 



= / p Tij ■=■ sin^-CW, - W a )t cos^ (W, + W 0 )t 

v I 

Using the assumptions on W ( and W 0 , the value of P 
Thus Equation 1.7 becomes 



5 = 




- X 2 - /2 



^ 27 zQ 

J SHR a - Ju 



dx + 



~V§UR (1 + y^JSR ) 



U27C £ 



- x*/2 



dx 



) 



- 00 



Finally, for ASK modulation, the optimum jammer is 

S(t) 



n.(t) -Jfr; 



IN 



and the probability of error becomes 



+ W 0 )T = 



, that Aj 
receiver 
and 1.7. 



is zero. 



given by 
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1 



where now 

P- 

and 



2_ ^ 

I27C 6 * 

(y/a* - /2JSR) 

2A, A „ 



-00 



- v/sNR (pi 

1 - x a /2 

v/27C 6 



+ V2JSR) 

dx 



.a. * z. 

A, + A e 



a-- 1 -p 



(A, - A 0 ) 

JL 
k. 0 



A, + A, 

The 'break point' for ASK occurs at JSR = 0^/2. Since qC < 
1, in terms of 'break point' efficiency, PSK is highest with 
a 'break point' occuring at a JSR of 1. FSK is next highest 
in efficiency with a 'break point’ occuring at a JSR of 1/2, 
and ASK is lowest in efficiency with a 'break point' 
occuring at JSR < 1/2. 



INCOHERENT RECEIVERS 



1. General 



the information 
exactly at the 



In the coherent systems considered, 
bearing signals were assumed to be known 
receiver . 

In noncoherent systems, the phase of the carrier 
signals is not available at the receiver so that the phase 
is treated as a random variable uniformly distributed over 
[0,27L]. As such, we may expect the performance of an 
incoherent system to be degraded in comparison to the 
performance of the corresponding coherent system. However, 
because of their simplicity, incoherent systems are widely 
used in many applications. 
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2. ASK ( On ^ Off key ing ) 

In the case of incoherent ASK (On - Off keying), the 
signals at the front-end of the receiver in the presence of 
noise only can be mathematicaly modeled by either 



H, : 


r(t) = A sin (W c t + 0) + n(t) 


0 4t 


or 

H 0 : 


r(t) = n(t) 


0 irt at 



The amplitude A, the frequency Wc and the time of arrival 
are assumed to be known except that the phase O is modeled 
as a random variable having an a priori density function 

0 ± 6 6. 2.7Z. 

otherwise 

The additive noise is again assumed to be zero mean, white 
Gaussian, with power spectral density level N 0 /2 watts/Hz. 

In the absence of a jamming waveform n.(t), the 
optimum receiver for decoding the binary information (which 
is transmitted via the use of either a sinusoid that is 
incoherently received, or no signal at all) in each interval 
of duration T sec, is the well-known quadrature receiver. 
Its structure is shown in Figure 6.2 [Ref. 2]. 

There is an alternate form of the quadrature 
receiver obtained by replacing the correlators with matched 
filters having an impulse response given by h(t) = sin Wc(T 
- t) and h(t) = cos Wc(T - t), 0 < t < T, and sampling the 
outputs at t = T (Figure 6.3) [Ref. 2]. 

Still another important alternate form of the quad- 
rature receiver is an incoherent matched filter followed by 
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an envelope detector and a sampler as shown in Figure 6.4 
[Ref. 2]. 

3. FSK 

For binary frequency shift keying (FSK), the 
received signals at the front end of the receiver in the 
presence of noise only can be mathematically modeled by 
either 

H , : r(t) = A sin (W,t + (0 ) + 0^.t^T 

or 

H 0 : r(t) = A sin (W 0 t + C|) ) + n ( t ) • o^t^T 

where the additive noise is again assumed to be zero mean, 
white Gaussian, with power spectral density level N 0 /2 
watts/Hz. Usually the frequencies W ( and W 0 differ substan- 
tially. The phases 0 and are statistically independent 
random variables, uniformly distributed over the interval 

[0,27t] . 

Insofar as the signal amplitudes is concerned, two 
cases are distinguished. In the first case, the information 
bearing signal amplitudes are known and equal. In the 
second case, the information bearing signal amplitudes are 
treated as independent random variables that are Rayleigh 
distributed due to multipath propagation. For this second 
case, the following hypothesis testing theoretic model is 
used 

H, : r(t ) = A sin (W,t + Q) + n(t) 0 



and 



H 0 : r(t) = B sin (W 0 t + (f) ) + n(t) O^t^T 
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where we assume that the signal amplitudes even though 
random, remain constant over a T - second interval. The 
noise is modeled as in the previous analysis, and the signal 
amplitudes have (a priori) probability density functions 
given by 

f A (a) = ex P 77S lzo 

and 

b , 2 - 

f B (b) = -p exp (- b -° 

where A| = E{ A 2 } = E{ B 2 } . 

The results associated with ASK incoherent 
receivers can be applied toward obtaining the optimum 
incoherent FSK receiver. This is worked out in Reference 2. 
The resulting receiver structure is shown in Figure 6.5 . 
Individual channel matched filters shown in Figure 6.5 can 
be essentially replaced with tuned bandpass filters having 
center frequencies corresponding to the "mark" and "space" 
frequencies W ( and W«, 



20 



II. ANALYSIS OF JAMMING ON ASK 



A. GENERAL 

In this chapter, the effect of a deterministic jammer 
waveform on the performance of an incoherent ASK receiver 
will be investigated. 

The signal at the front-end of the receiver is given by 

Equation 1.8 with the modification that under either 

hypothesis, a jamming waveform n ; (t) is present during the 

o 

time interval [0,T]. 

For the coherent ASK receiver, the optimum energy 
constrained jammer waveform could be obtained by using 
Equation 1.5 . However, the complexity of the mathematical 
expression for probability of error in the noncoherent case 
makes it very difficult if not impossible to derive the 
optimum energy constrained jammer waveform in closed form as 
will be seen in the sequel. 

A reasonable postulation is that the optimum jamming 
waveform for the coherent ASK receiver can act as a good 
jammer for the incoherent ASK receiver also. Thus, such 
near optimum jammer signals are studied and evaluated in 
terms of their effect on the performance of the noncoherent 
ASK receiver. 



B. ANALYSIS WITH NEAR OPTIMUM JAMMER 

Analysis of the incoherent receiver starts from the 
mathematical model of the receiver front-end input waveform 
r(t) given by either 

H, : r(t) = A sin (W c t + 0 ) + n(t) + n. (t) 0 £T 



or 
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0 it iT 



H a : r(t) = n(t) + n- (t) 

J 

where n.(t) is the jammer waveform present during the time 

v] 

interval [0 , T] . 

In the absence of n.(t), the optimum receiver for the 

binary ASK problem is well-known. Its derivation is well 

documented in the statistical detection theory literature 

[Ref. 2]. The receiver structure is shown in Figure 6.2 

In this section, the effect of n.(t) on this receiver is 

o 

analyzed by evaluating the resulting Pe , under the assump- 
tion that n.(t) is a deterministic jammer waveform, however 
unknown to the receiver itself. Receiver performance 
requires determination of the statistics of either G 2 or G, 
where G 2 is the output of the quadrature detector and is 
given by 

G = X + Y 

where 

X =£r(t) sin W c t dt = (r,S) 

Y = j T r(t) cos W c t dt = (r,C) 

Provided that the random variable Q) is fixed to some value 
9, X and Y are conditional Gaussian random variables with 

E | X ! H , ,0j = (AS e ,S) + (n^ ,S) = m >j9 

and 

E ’Qj = (AS 0 ,C) + (n. ,C) = m y|e 
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where S & is used in place of sin(Wc t + 0 ). It can also 
be seen that 

Var jx|H,,9j = E {(n,S) a } 



M,T 

4 



sin 2W C T 
2W t T • 



and 



= { ^|h 0 | 

s imil arly 




N 0 T 

4 



i + 



sin 2W^ T 
~2W C T 



=Var 




For convenience we assume WcT = n 7t , where n is an integer. 
Thus the sine terms above vanish resulting in var{X I (-{, > 0 } 
= var{Y | bli , & ) = N 0 T/4 = (J 1 . If this assumption is not 
made, an additional term results. However, if 7C /Wc << T, 
the additional sine term is small and may be neglected. 
Furthermore the covariance 



E 




E {x|H;.e}] • [ 

{(n,S) (n,C)J 



Y 



- e{y|h. ,e} ] 




M 0 T 

"V 



1 - cos 2W C T 
2W t T 



0 



i = 0,1 



provided the assumptions on Wc hold. This implies that for 
any given value of 0 , both X and Y are uncorrelated 
Gaussian random variables and therefore statistically inde- 
pendent. The density function of G 2 conditioned on the 
phase 0 is noncentral Chi-Squared distributed, and is given 
by 
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f e r i l0 ( §l H ' ’ © } = 37^ ex P^- Ic (^ 



2 (J 



, g >0 



and zero otherwise. The variance 2 is defined above and 
becomes 



S = + mrie = E Z {x|H,,9} + E* {y|H, , © | 

= ^(AS# ,S) + (n.,S)]*+ [(AS 0 ,C) + (n.,C)J* 



Due to our assumption WcT = n 7C , we obtain 

x 

S = \ ~ 2 ~) + AT [( n j>S) cos ® + ^ n j sin ©j 

+ (n ,S^ + (n ,C) a 



( 2 . 1 ) 



Similarly, if r(t) consists of noise and jammer only, then X 
and Y are also independent Gaussian random variables with 



E { X|H 0 f = (n. ,S) 



<4 

= m, 



and 



E | Y|H 0 J = (n J5 C) = m Y 

The density function of G 2 assuming no signal is sent is 

V «*(- * 7 $) 



2 (J 

and zero otherwise, where 

S = < + m y = (n. ,S) Z + ( n . ,cf 



(2.2) 
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The resulting average probability of error is given by 






Pp ' P(H °i + P " P(H ' } 71 1 

p ( H „) j f^Cg|H 0 ) dg + P(H, ) | f^a. 



Cg|H, ) dg 



where Tj is the threshold with which G is compared in order 
to decide which is the true hypothesis. Observe that 



V ( g|H,) = Vie ( Sl H ” © > 

The first integral in the above expression for Pe is 
(probability of false alarm) . It can be expressed as 
follows 



•CO 

J_ g + s 

2o *exp (- 



d g 



T 



-CO 

V exp (- 1 1 21 ) I ( q'v ) dv 



'la 

- Q ( a' , n /a ) 



where 

r 00 i a. 

Q ( a , B ) = V exp (- — * - ) 

J /3 

is a well-known tabulated function 
Q- function. The second integral in the 
Pe is P M (probability of miss). It 
follows 



I 0 ( a V) dV 

called the Marcum 
above expression for 
can be expressed as 




H, ,9) f 0 ( 0 ) d 0 



dg 



( 2 . 3 ) 
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fco 




C-J f 1 


— 


f (0) 




If (g|H, ,0) dg 




6 




G | 0 


> 


■CO 


_ J 


0 



where in the second equality the order of the integration 
has been changed. The inner integral of Equation 2.3 can be 
expressed in term of the Marcum Q-function as 



L_L_ exp ( - g ) I ({El.) dg 

2o J 2 a 1 a 2 



fV/a 



v exp (- v + a ) I 0 c ctV) dV 



= 1 - Q ( a , n / a ) 



where q( = J~S "/ (j . Therefore the probability of error can be 



written as 



/s 7 



P e = P(H e ) Q n /a ) 



(2.4) 



+ P(H, ) 



(2. 7Z 



1 - 



2 IT 



Js 



Q ( , n /a )* d 9 






where the dependence on 0 is imbedded in the term / S and 
the threshold Tj for the jammer absent case is obtained as 
the solution to the equation 



- A T/ 2 N 0 

g I 0 (, 2 An/N 0 ) = r 



P(H 0 ) 
’ R " PCH, ) 



which can be equivalently put in the form 



- A 2 " T/ 2 N, 



e 



i 0 [Jir™. ( n / 



a) ] 



= R 



The I 0 ( • ) function used here and also previously used in 
conjunction with the development of Pe is the modified 
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Bessel function of the first kind. The appropriate setting 
of the threshold is an important issue for both coherent as 
well as noncoherent receivers. One could use the threshold 
setting that would be derived from the analysis of receivers 
operating in additive white Gaussian noise interference 
only. This approach can be quite unsatisfactory as demon- 
strated in Reference 1. A better approach would be to 
obtain Pe as a function of the threshold, and then search 
for the threshold that minimizes Pe. While this approach is 
intuitively appealing, it is often mathematically intrac- 
table. While the threshold issue is not addressed in this 
particular section, it will be discussed in more detail in 
the sequel and simulation results will be presented. 

If the definition of average signal energy previously 
introduced is used, we have E = A 2 T/4 which is reduced in 
half in comparison to the binary signal transmission case, 
due to the fact that the information bearing signal(s) do 
not have equal energy. In order to afford comparisons with 
the coherent receiver case, we will implicitely boost the 
value of signal amplitude A, to obtain E = A 2 T/2 in order to 
have agreement with previous cases insofar as signal ener- 
gies is concerned. Thus the threshold determination equation 
now becomes 




-SNR 



( n / a ) ] = R 



P(H, ) 
= F(H, ) 



(2.5) 



If we assume P(H,) = P(H„) = 1/2, then 




de + Q 




and the threshold setting equation becomes 



-SNR 




e 



1 
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If we are to find the optimum jammer waveform so as to 
maximize Pe , an attempt must be made to solve 

bp ap 

— = 0 and as-= o 

Unfortunately the resultant equations are mathematically 
involved and do not appear readily solvable for n.(t). It 

sj 

seems however that a 'good' jammer waveform can be postu- 
lated based on the results obtained for coherent ASK. It was 

found for that case that the optimum n.(t) is a tone at the 

o 

carrier frequency. Thus the following jammer waveform can be 
used as a potential near optimum jammer, namely 

f~p ( 2 6 ) 

n J^ = Y 1 1’J'r s; ’- n W cP , 0£tiT 

Observe that with this choice, | | n - | | 2 = P-*. . The prob- 

J 0 

ability of error Pe can now be determined using the 
threshold setting equation (Eqn. 2.5) and the previously 
derived expressions for S and S'. The effect of the near 
optimum jammer waveform on the receiver (i.e., incoherent 
receiver performance) can be analyzed by evaluating Pe as a 
function of n.(t) using Equation 2.6 . Note that in 

si 

Equation 2.6 the jammer energy is P,,. . It can be shown that 

J 

<n J> S) = w «t 

(nj,C) =J T sin W c t cos W c t dt = 0 



Then from Equations 2.1 and 2.2, the probability of error Pe 
given by Equation 2.4 can be expressed in terms of SNR 
(signal- to-noise ratio) and JSR (j ammer- to- signal ratio) 
using the fact that 
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~dj = 2 SNR (1 + 2 /jSR cos© + JSR) 

= 2 SNR- JSR 

and the fact that Tj j (j can be obtained from Equation 2.5. 

For equally likely hypotheses, the probability of 
receiver error (Equation 2.4) can be expressed by 



Pp = ^-[l +^Q ( 72SNR JSR , ^/(J) 

i- /"o^ ( / 2 SNR (1 + 2,/JSR cosQ + JSR ) 



7 ]/(j'> d 0] 



(2.7) 



From Equation 2.7, Pe greater than 0.5 can occur if 



Q (ct ,3 ) > 2 “ Q U',B > d6 



( 2 . 8 ) 



where 



d'- 72SNR (1 + 2/lSR cos0+ JSR) 

Ct - J 2SNR~JSR and ^ /(j 

Since the periodic sinusoidal function imbedded in 7 ' is 
integrated for one period, the condition p( ^ ' might 

cause the condition of Equation 2.8 to be satisfied for a 
fixed value of ^ . Thus there is a possibility of obtaining 
Pe greater than 0.5 for a JSR value beyond the critical 
value of 0.25 which can be obtained from the condition 
oi ~ oL ’ > that is , 

/2SHR-JSR = v/2SNR (1 + 2~/jsr cos0 + JSR) (2.9) 

7 JSR cosQ = — -L~ 



29 



Receiver performance in the presence of a jammer is 
expressed in terms of signal- to-noise ratio (SNR) and 
j ammer- to- s iganl ratio (JSR) . A rectangular pulse of dura- 
tion T seconds has amplitude spectrum AT sine Tf, and B = 
2/T is a rough measure of its bandwidth. Thus the expression 
of the form 



A 2 T/2 



A Z /2 



A z /z 



( 2 . 10 ) 



Jh 2 PSD-B 
2 T 



can be interpreted as the ratio of signal power to noise 
power (SNR) in the signal bandwidth [Ref. 5]. The term JSR 
can be expressed by 



_ jammer energy 



( 2 . 11 ) 



A i T/2 signal energy 



where P 71 . represents the jammer energy defined before 

■J 



C. VARIABLE THRESHOLDING EFFECT 



It is apparent that when a fixed threshold value is used 
by the receiver, the effect of the jammer waveform having 
energy such that the resulting JSR value is above a certain 
level, is such that the receiver may be rendered inoperable. 

It could therefore be suggested that setting the 
threshold based on the value P(H|) or P(H 0 ) only (see Eqn. 
2.5), may not be desireable. Values of threshold other than 
some fixed value may result in improved receiver perform- 
ance, in other words, reduce the jamming effect. 

Recall that the threshold setting equation is of the 

form 



-SNR 

0 I G ( y 2SNR a) = R 



( 2 . 12 ) 
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where o( - 71 / (J . we now can attempt to reduce Pe by an 
appropriate choice of threshold ^ . By varying the value of 
R (instead of using R = 1) we can obtain this threshold 
which denote ' from the expression 

I SNR 

a' = — r — i;'(r e > 

v/2SNR 

where I 0 ' l (X) is the inverse modified Bessel function of the 
first kind. Using an approximation to I 0 (X) which is given 
by 

P x 

loCx) = ~ for x^l , 

Equation 2.12 can be rewritten as 



. (2.13) 

4 y 2 SNR a' - 2 

- 2 2tt + J^rl( 2 SNR) + 4 SNR + i'N R 

From Equation 2.13 , it can be recognized that when SNR is 
large, variation of the value of R does not significantly 
affect the threshold value qC ' because taking the logarithm 
of R reduces the effect of R further. Thus the large value 
of SNR suppresses the effect of variation of the term J2n R 
on ^ 1 . This limited variable thresholding effect on the 
receiver which is under significant jamming is analyzed in 
graphical form in Chapter 5 for various values of R. 
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III. ANALYSIS OF JAMMING ON FSK 



A. GENERAL 

For binary incoherent FSK with a jammer present, the 
received signals under the two hypotheses are either 

H, : r(t) = A sin (W,t + 0 ) + n(t) + n.(t) 0 

j 

or 

H 0 : r(t) = A sin (W 0 t + (j)) + n(t) + nj(t) 0 

By separating the frequencies W, and W 0 sufficiently, we can 
form signals that are orthogonal, have equal energy, and 
have the same advantage of ease of generation. 

The modified FSK receiver structure which is capable of 
varying the output of the each envelope detector (which is 
followed by a multiplier) is diagramed in Figure 6.6 . The 
optimum receiver for the case where no jammer is present can 
be derived from statistical decision theory and is shown in 
Figure 6.5 . This receiver can be obtained by setting = 

1/2 in the modified FSK receiver shown in Figure 6.6 . 

In practice, incoherent FSK is widely used because of 
its simple receiver structure, its small performance penalty 
due to lack of phase coherence and its more efficient use of 
signal energies in comparison to incoherent ASK. In addi- 
tion, we are not faced with the difficulty imposed by a 
threshold that must change with SNR as is the case with 
incoherent ASK. Thus, a receiver which is known to be 
optimum for incoherent FSK transmission, has been modified 
by including some channel weighting. This has been done in 
order to be able to determine whether or not such channel 
weighting can reduce the effect of the jammer. 
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This chapter is devoted to investigating the performance 
of the modified incoherent FSK receiver in the presence of 
jamming and additive white Gaussian noise. By letting ~ 
1/2 we can as a byproduct obtain the performance of the 
conventional incoherent FSK receiver of Figure 6.5 in the 
presence of jamming and additive white Gaussian noise. 



B. ANALYSIS WITH NEAR OPTIMUM JAMMER 

The modified receiver performance can be obtained by 
introducing a null hypothesis (no signal) as a dummy 
hypothesis and by following the same reasoning as in the 
analysis of incoherent ASK presented in the previous 
chapter . 

The receiver function is to compare the envelopes at the 
output of each channel once every T seconds and decide in 
favor of the larger of the two envelopes (Figure 6.6). For 
the purpose of analysis, let us first assume that a 'mark' 
signal has been transmitted, that is, the hypothesis H, is 
assumed to be true. An error is committed if V 0 exceeds V| . 
An error is also committed if V t is larger than V $ when a 
'space' signal has been transmitted, that is, the hypothesis 
H 0 is assumed to be true [Ref. 3]. 

Let Pe, denote the probability of the first type of 
error described above, which is expressed as Prob.(V,;, > V ( | 
Hj ) . Under the assumption that a 'mark' signal has been 
sent, the output q , of one of the envelope detectors is 
given by 



qf = + Y, 2 - 



where 



X, = r(t) sin W,t dt ; Y, = / r(t) cos W.t dt 
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Observe that Xj and Y ( conditioned on the phase and either 
of the two hypotheses are Gaussian random variables with 



E {x,| H , , 0j = JT A sin (W,t + 0 ) sin W,t dt 

+ f T n .(t) sin W,t dt = (S,,S), +(n.,S), 

J o < J 



and 



E 



cos W,t 



{y, | H, ,0 j = j A sin (W, t + Q ) 

fr Jo 

+ I n.(t) cos W, t dt = (S, ,C), + 

Jo J 



dt 

(n.,C), 

J 



where Sj represents the function A sin ( W , t + 0 ), S repre- 
sents the function sin W| t and C represents the function cos 
W ( t. Likewise, assuming again that W ( T = n7L, n an integer, 
and that n(t) is zero mean white Gaussian noise with PSD 
level No/2 watts/Hz, we obtain 

Var Jx,|H, ,©j = Var |y,|H ( ,0J = 

Furthermore, it can be shown that 

E |[x, - E{x,|H,,e}](Y l - E{Y,|H,,e}]|H,,eJ = 0 

so that the conditional r.v.'s X, and Y ( are uncorrelated 
and therefore independent. The sum involving random vari- 
ables X, and Y|, and producing q 2 , will result in a non- 
central Chi-Squared distribution so that 



f oti0 Cq,|H,,e) = 2 a A exp (-V;dr ) In 



2 a- 



■° 0 2 - 



q, 



where 



S„ = E Z {x,|H, ,©j + E* { Y, | H, ,©J 

= [ ( S, ,S), + (n.,S),]* + [CS, ,C), + (n^C),] 
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and 0 * 1 - N 0 T/4. Using standard random variable transforma- 

tion techniques, it is not difficult to show that 



so 



that 




q* + s, 

2 o 1 



) i, 




a 



q, ^ o 



(3.1) 

( q, I H, ) = f aAd (q, [ H, ,0) (0) d0 

= 1 PS exp/ qf + S„ \ i o ( q, /sT ) d p> 

2tt|~1 \ 2 a*- J a 1 

-/o O 



where the dependence on 0 is imbedded in the term S l( . We 
now need to obtain the statistics of the output of the 
multiplier following the upper envelope detector. (See 
Figure 6 . 6 ). That is, the probability density function 
f v (V, | H,) needs to be derived. From standard transforma- 
tion theory, using the relation V, = 2(1 - )q ( , it can be 
shown that 



f Vi (V,|H, ) = 






v, 



2(1- a) 2 ( 1- a ) 

dzTC 

Vi 



|h.) 



(3.2) 



4 7t ( 1- a ) 
I 



2 ( 1- a ) exp ( 



( 2 ( 1 - a )) * S " ' , 



0 a 

(TTiVV^ ) 

' 77X ' 



2 o' 



d8 



On the other hand, the output of the lower envelope 
detector when H | is assumed to be the true hypothesis, is 
given by 
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o 



= xS Y * 
o 0 



where 



X 0 = r(t) sin W.t dt 



Y. 



r(t) cos VLt dt 



Following a similiar procedure as used above, it can be 
shown that 



e|y o |H, , 8 j = j A sin (W, t + 0 ) cos W c t dt 

Jo 



+ ln.(t) cos Q W t dt = (S , ,C) + (n. ,C) 



and 



t + 0 ) sin W 0 1 dt 



, 8 1 = j A sin (W, 

so 

+ J nj (t) sin 0 W t dt = (S, ,S) C + (n. ,S) a 

s n 



Here, S and C have the same meaning previously defined 
except that the subscript "0" outside the inner products 
implies that we should interpret S as sin W c t and C as cos 
W 0 t . It can also be demonstrated that 

Var | X 0 |H,,q} = Van |y 0 |H ( ,©} =-^L 
and also that 



E [ [ X 0 - E [ XjH,,©} ] [ Y 0 - E ( Y 0 |H,,eJ]|H,,ej= 0 

so that the conditional r.v.'s X 0 and Y 0 are uncorrelated, 
hence independent. Thus similiar to Equation 3.1, the 
expression for the conditional density function of q 0 
becomes 
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- 2 *- 

<lo 



Q 0 (q ^ H,) " 2-n I aZex P ( 



q^+ S,. 

2o z 






d 8 



q 0 ^° 



where 



'0\ 



= E a { x o |h m 0J +e 2 {y 0 |h m gJ 

= [ CS, ,S) 0 + (rij,S) £) ] i + [(S,,c ) 0 + Cn.,C) fl ] 



and dependence on 6 now is imbedded in the term S 0) . 
Applying again the random variable transformation V D = 
2^q 0 , the conditional density function of Vo is given by 



f Vo (V 0 |H,) 





(3.3) 



Vo _ 

( 2 a7 S °‘ ^d8 



In order to compute the probability of error, we can now 
use the previous expressions for the conditional probability 
density functions which are derived assuming a 'mark' signal 
has been sent. That is, for a given value of V| , an error 
is made if V 0 > V|. Thus the average error probability is 
found by averaging the conditional error probability given 
by 



r* 

P(V, ) = | f (V 0 1 H, ,V, ) dV a 



over all Vi . That is, 



P e , = Prob. £ V 0 > V,| H , j 



(3.4) 
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(Vo | H , ) dV 0 



f Vi (V,|h, ) dV, 



Substituting for the density functions in Equation 3.4, from 
Equations 3.2 and 3.3, we obtain 



P 

e. 






2 7C ) 0 L 




° 

,, ° aeX P (" 




(3.5) 



Io 




fy( V | | H , ) dV, 



where the order of the integration has been changed for 
computational convenience. In the above equation the inner- 
most integral can be expressed in terms of the Marcum 
Q- function as follows 




exp 



(- 




I e ( 2a° dVo 
o^ 



\fsZ ^ V, 
a 



■ Q < ^ > 



Then, Equation 3.5 for Pe becomes 



ei 



1 

2 TT -'0 



ui 



! q < ^ de 



f Vi (V.|H,) dV, 



(3.6) 



(iX. 



Q ( 



\fs^ 



2 TT Jo 



1-q 

a 



x ) de 



Cx n. 



I x exp ( - 

2 TT Jo 



X 

X + 



S ,, 



) i e ( x de 



dx 
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where in the second equality the change of variable x = 

V , has been m3de . 

2 ( l-ot)cr 

From the orthogonality property of the signal pair used 
(which is obtained by assuming sufficient separation between 
two frequencies and that W ( as well as W Q are large), we 
have 

j sin (W,t + 0 ) sin W.t dt = 0 

J o 

so that the term S 0 | is independent of Q . Therefore 
Equation 3.6 can be rewritten in the following form 



P = 
e, 




/F°l s x) x exp ( - 

o ’ a 



x* + 




2 . 




2 



(3.7) 



I 0 (x 



•J 3, « 



) dx 



d9 



where the order of integration has been changed. 
Furthermore, using the following formula involving an inte- 
gral of a Marcum Q-function [Ref. A] 



f . o.i Ri . R, / cti ^ R| \ j / ctiR, > 

h ( 157’ -57 * 557 exp “W - ' ^ 

Jq 



dR, 



o,+ V 



1 - Q 



a, 



*+ b* 



o" + a, » V a i 



X 

o% 



G, i + aj- 



Q < 



2. 

2. 



X 

a, 



T a a ’ vdf + of- 



fa, 3 - + 



the inner integral in Equation 3.7 can be simplified in such 
a way that Pe, becomes 

( 3 . 8 ) 



p = 



(1- a) 



e ' (1- of +ct* 



r i 


rxJt 


yin" j d-d ) y"so( / ci x ^ 


( 

i 

4 


'w' 

O' 

o 


u V( 1- a f + of , a \/(l-o . 
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+ 



z 



a 

( 1-a f+a* 



‘ 1 

2 TT 



p-7 C 

Q ( 



v^Soi /- a 1 /Si7 / (1-a ) . , Q 

a '/Cl-a ) +a 1 , “V(l-a r+a 2 - ) d9 



where the dependence on 0 is imbedded in the term S (l . 

Following exactly the same procedure used in obtaining 
the expression for Pe ( , it can be established that the 
expression for Pe 0 , which denotes the error probability when 
H 0 is assumed , that is , Prob . ( V t > V 0 | H 0 ) takes the form 



(3.9) 



P = 

6 o 



a 



(1- oO x +a x 



1 r sz / n- qf 
7F J 0 ^ ( o /l-aT+ a 1 



v /ST / (l-cQ x 
(T V 




(1- a) 
(1- af+a 1 



1 
2 TT 



I n , \/s<»* / ( 1 — CL ) ^/Soo / a " ” \ J /f 

I Q ( — '/(T-a )W /(“ ^ Aa 1 



where only the term S oo is dependent on . Therefore the 
total average probability of error (Pe) can be obtained from 



p e = -g- ( P e , + P eo ) 



assuming that the two hypotheses are equally likely. 

For the special case ^ = 1/2, the performance of the 

optimum receiver in the presence of jamming is given by 

i - 4- Hq < ^ a- & 



P e Cc =A) 



2 TT 



'CT~ o , 'TT a 



•) 



(3.10) 



- Q ( vt IF,/? IF’} oe 

_1_ fi* f 

2 TT / / Q ( — 1 '/S 7 ff 

A r ^ a , rr 



_i_ ys, 0 1 /s,J . ] 
■2 o ,/2 —f^l 
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~ Q ( 



Observe that Equation 3.10 yields the performance of a 
conventional incoherent FSK receiver (with no channel 
weighting) in the presence of jamming. If we now use as the 
jamming waveform the jammer which is optimum against a cohe- 
rent FSK receiver, that is 



n .(t) 
J 



= sin \ <W, - H)t cos j(W, + 

--IK 

J T L 



W fl )t 



sin W,t - sin 



W 0 t 



( 3 . 11 ) 



then the terms S ;k (i,k=0,l) in Equations 3.8 and 3.9, which 
are a fund 
as follows 



are a function of the jammer waveform n.(t), can be computed 

0 



and 



a 5 

= £ (S , ,S), + ( n j , S ) , + £(S l ,C), + (nj ,c) t J 



♦ _A jJM cos8 ♦ V 1 

s 2 11 



2 

S oi = [ (S, ,S) Q + (nj,C) 0 J + [ ( s, ,C ) 0 + (n. ,c)J 



■P.jT 



= ( 



[(S P ,S ) 0 + (n, ,C) 0 J + [(S 0 ,C ) 0 + (n.,C) 0 _ 

4Lf. cos 4 + 

* 2 4 - 



j(S 0 ,S), + (n.,S)J + [<S 0 ,C) ( + (n.,C)] 



4 - 
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Thus the probability of error of the unmodified receiver 
(Equation 3.10) can be expressed in terms of SNR and JSR 
only (defined by Equations 2.10 and 2.11) as follows 







where 



(/.»= SNR (2 + 2 /2JSR' cos 0 +JSP) 
oloo- SNR (2-2 / 2 JSR cos (f> + JSF) 
f?„,= fio = SNR JSR 



Receiver performance can now be evaluated as a function of 
SNR for fixed values of JSR. 

To provide further insight into the performance of the 
modified incoherent receiver, the effect of varying can 
be analyzed via computation of the probability of error Pe 
from Equations 3.8 and 3.9. In terms of SNR and JSR, Pe 
becomes 

P e = P(H , ) P e , + P(H 0 ) P eo (3 13) 

= P(H.) C,[ 1 - J'q' de] • 

+ P(H|) de] 

+ P(H 0 ) C^l - 2 ^ f Q (/Ci c 77o ,Jc t @, 9 ) d<J> 

+ P(H 0 ) C,[ — jr J ^ Q (J c p ig o4t>) dcf 
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where 



P . ( 1 - r = ** 

' U-oifroC 3, ’ Z O-eON-^ 3, 

The effect of a jammer that uses a single tone only to jam 
either 'mark' or 'space' signals, can be analyzed by setting 



n. (t) 

j 

or 




sin W,t 



n. (t) 

J 




(3.14) 



If the first choice is used, we obtain 



(3.15) 



Ig = P(H,) C, [ 1 - 27 hJ[q ) d9 ] 

+ P(H, ) C 2 [ 0»,' ,/cTZT ) ds] 

+ P(H 0 ) C 2 [ 1 - C/fTTTT ,JZJ^ ) d<|>] 

+ P(H ^ Ci [ (v/ ^ r H 



where 

o/„= 2 SNR (1 + 2/jSR cos 9 + JSR) 
(%l'= o 
c/.*= 2 SNR 
ft/= 2 SNR JSR 
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and if the second choice is used, we obtain 



p e = p(h , ) c, [ i - J-$ (Jc, Jc>p 'o ) de] 
+ p(h, ) c 2 [ jjf Ja ) de] 

+ P(H e ) C z [ 1 - -±-\q Jc7p>', ) dcp] 

+ P(H 0 ) C, [ Q ,/cT7 T ) dcp] 



The effect of varying the threshold oC on receiver perform- 
ance for the jammer waveform of Equation 3.14 is discussed 
in Chapter 5. 



C. FREQUENCY MODULATION SWEEP JAMMING 



In certain situations, the need to jam a certain 
frequency band rather than a discrete set of frequencies 
using tone jammers may arise. Therefore an FM sweep jammer 
will be proposed, analyzed, and its effect on the incoherent 
FSK receiver investigated as a function of the number of 
times the jammer sweeps the signal band during the signaling 
interval [0,T]. The mathematical model used for an FM sweep 
jammer is 

n j ( t ) = sin [ W s t + K 4 |a cos VLt dt ] 0^t£T 
After carring out the integration we obtain 



n. (t ) 

j 




sin 



W s t 






sin 



W.t 

J 



+ 9 
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where p = K_p&/Wj and 0 is a deterministic phase angle. 
The instantaneous jammer frequency is 



W i ( t) = W s + R w. cos W.t 
r J J 

so that the instantaneous jammer frequency covers the band 



Assume that WsT = 2 7C 1 and W. T = 
3 o 

In order to determine 



(W 5 - (3Wj , W s ♦ 0Wj). 

27 Ck, where 1 and k are integers, 
receiver performance in the presence of such a jammer, the 
parameters S; K (i,k = 0,1) which are a function of n.(t), 
must be determined. Before so doing, we evaluate 



(nj.S)* 




sin ( Wgt + |g sin W.t) sin W k t dt (3 16) 



K = 0 , 1 

where since the fixed phase 0 represents a time delay only, 
it is set to zero for computational ease. Equation 3.16 can 
be rewritten as follows 




(3.17) 



By using the well-known Bessel function coefficient expan- 
sion for each cosine term in Equation 3.17, the following is 
obtained 



H 

. S?J(Ws 


r SirJ 


[(tu s -a> K )+nia]T 


(w s -u/ k ) + 

+u/0 + / w.1t ' 


iUV- 

J 


(Ws 


+ UUk) t H LU. 


- 



(3.18) 
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Since the FSK signal covers approximately the band (W 0 " 
4 7C/T , W, + 47C/T) (assuming W t > W 9 ), it is reasonable to 
set 

W 5 = -i-( W, + w 0 ) (3.19) 



so that the instantaneous jammer frequency band (W$ - (3 Wj , 
W s + @ Wj ) completely covers the signal band provided that 



W< 



- (3W_. = w 0 - -A* 



W, 



+ £ W J 



= V? 



4-7C 



This means that 

pw. = -^-( W, - W„) + 

must be satisfied. Since for FSK signaling we have assumed 
that 



( W, - w 0 ) T = 17C ; (W, + w e ) T = m7L 

where 1 and m are integers, and W*T = 2Xk, it is apparent 

J 

that 

r- tH- + 1/2] 

Note now that the integer k determines the number of times 
the jamming waveform will sweep the signal band in one bit 
interval. Thus from Equation 3.18 and Equation 3.19, we 
have 

(n. ,S)i 
J 



which becomes 



SI ? j„< & )i 


\ StT) 


(i(w.- w,)+7HV.]T 
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(3.20) 



P T 2 

n. I 



(n. ,S), 
J 



I J « ( f 3 



r-.-co 



f Sfn nfenk-f/?) 

7^(2^ K-J?/a) 



Sin ±l/Z} 

7i^2nk-nn + j2/£) 



and 



(n.,S) e 

j 



fis 



£T 



IJ.' P 



T's -OP 



s j n t3 1^ 0 ) j- t 

+3^4-711^. 



SinfxC^ 1 ~ lip 4-71^.] T 

2, ( — W-A) UJ* 



which becomes 



(n. ,S ) 0 
o 



2k 



rT 



p I J - ( P 



SiTi 7 r( 2 rm-hl/£) 
7L(anK -t ^/2) 
_ Sin 7C.(anK-f7H-^) 

7i^ar\K + ^ -^/a) 



(3.21) 



If 1 is an even integer, then Equation 3.20 and Equation 
3.21 at most contain two terms respectively. These terms can 
exist only if the argument of all SINC functions is zero. 
Therefore Equation 3.20 and Equation 3.21 can be simplified 
to yield 

<n,,S), =/fg£ [j„,< (3 > - < fn] 



where n, = l/4k, = - (m + 1/2) / 2k and 

(n.,S). =/E. [j ( (3 ) - J < (3 )] 

J t— ~ 

where n^ = - l/4k, n^ = - (m - 1/2) / 2k. 



In view of practical communication system constraints, 
the integer m will typically be much larger than 1 because m 
and 1 represent the sum and difference of the signaling 
frequencies respectively. In particular we see that 
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J K ( ft )<< 1 for j<3 << k. Hence the terms J r1x ( p ) and 

( A ) can be neglected so that 



and 




(3.22) 



(3.23) 



The other parameters to be considered are 

(n.,C) K = ( fivSsin (W 5 t + Q sin W.t) cos W K t dt 

J V ^ T r j 

where |3 = K / Wj . Using again the Bessel function coef- 
ficient expansion and appropriate trigonometric identities, 
we obtain 



(n.,C), 

j 1 



•m-T 



I f ^ 

m - -vn 



-g(2nK+in + X/z'j 



s;Vt(2TiK-V2) 

-f(ZrlK- X/2) 



(3.24) 



and 



/ P - T ^ 

(n , C) 0 =y-~ ^ J„ ( 0 ) 



W- -CO 



sin i -?f'( 2nk + 1/2) 
( 2nk + 1/2) 



(3.25) 



sin 4 -J-(2nk + m -1/2) 
"( 2nk + m -1/2 ) 
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We note that all the terms in Equation 3 . 24 and Equation 
3.25 take on the form of a SINC function times a sine func- 
tion. We can readily show that for even values of m and 1, 



( n J ’ C ) , = ( n . , C) 0 = 0 



(3.26) 



Using the mathematical forms given by Equations 3.22, 3.23 
and 3.26, the performance of the receiver with FH sweep 
jamming is obtained via the use of Equations 3.8 and 3.9 and 
evaluated from the expression 



(3.27) 

P. = P(H,) c, [l - c/ZT^ ,/ZTpZ > de] 

+ PCH,) c 2 Jj 1 (Jc x p ol jtz ) de] 

+ PCHo) c * l 1 ~ ) d<p] 

+ P(Ho) C| [~Jtl ( /^ > dtp] 



where 



<4= 2 SNR [ 1 ♦ 2m J „, ( (3 ) cos0 ♦ JSR ( p ,J 

c4.= 2 StlR [ 1 . 2 /35 r j, 3 ( p , cos<( , t JSR ^ ( p j 

(3>i = (3,0 = 2 SMR JSR J^ ( p ) 
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C t and Cg. have been defined in Equation 3.13. As before, in 
this equation SNR and JSR are defined by Equations 2.10 and 
2.11 respectively. 

The performance results to be presented are a function 
of SNR, JSR, and the number of times the jammer sweeps the 
signal band. This is discussed in greater detail in Chapter 
5. 
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IV. ANALYSIS OF FSK IN THE PRESENCE OF JAMMING AND FADING 



A. GENERAL 

In certain propagation media, the received signal trans- 
mitted via a free space channel is often subject to fading, 
a phenomenon caused by multipath propagation and the equiva- 
lent addition of random phasors. The vector sum signal will 
have an envelope which changes with time resulting in an 
effect known as fading. 

The fading signal model often utilized is that of the 
nonselective , slow fading, Rayleigh distributed signal 
amplitude where it is assumed that the amplitude, while 
random, remains constant over the time interval [0,T]. Thus 
the received envelope is now random with probability density 
function 

0 2 - 

f A (a) = — |-exp(-4-) a > 0 

A E E 

where E = E{A 2 } denotes the mean squared value of the signal 
amplitude . 

The optimum (minimum probability of error) receiver 
structure in the presence of additive white Gaussian noise 
(however with no jamming present) is the same as that for 
noncoherent nonfading FSK and its structure is therefore 
given by the receiver of Figure 5. 

B. ANALYSIS WITH NEAR OPTIMUM JAMMER 

We shall analyze the effect of Rayleigh fading on an 
incoherent FSK receiver operating in the presence of a 
jammer in addition to the additive white Gaussian noise. 
Thus, the signals at the front end of the receiver are 
either 
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or 



H, : r(t) = A sin (W, t +0) + n(t) + n.(t) O^t^T 

J “ 



H : r(t) = B sin (V7 0 1 + (J) ) + n(t) + n.(t) 0^t^.T 

J 



where the amplitudes A and B are independent identically 
distributed random variables having density functions 



f A (a) = i^ exp(- 2A 1 ) a£0 

f B (t.) = -^laxp(- bZO 



with E{ A) = E{B } =J -g A 0 and E{A 2 > = E(B 2 ) = 2A 2 . 

Note that FSK with fading is the same as what was previ- 
ously analyzed (FSK without fading) except that now the 
signal amplitudes are random varibles. Therefore for a 
fixed value of A, the error probability when a "mark" signal 
has been transmitted is the same as that of the noncoherent 
(nonfading) FSK case. That is, we can use the results of 
the preceding chapter and in particular make use of Equation 
3.7, to obtain (for the special case of ^ = 1/2) 



P (A) 

e, 





x - + fHL\ 

x exp C - v a / ) 
2 



I 0 (x — ) dx 



d 8 



where the dependence on A is due to the fact that the term 
S n defined previously depends on A. Then, P e becomes 

foo * 
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exp ( - 



:a; 



:) dA 



J*C x 1 

Q(^d-,x) x exp( — j) 



& 



ysj exp 



( -forrf-A* n 



2 A 



I 0 (xf, dA | dx] 



d0 



where the last equality has been obtained by interchanging 
the order of integration. This equation does not appear to 
be readily simplif iable . If we take a closer look at the 



innermost integral in Equation 4.1, 
evaluated using 



tU 1 I y (ct t) exp(-p 2 t 2 ) dt 



that integral might be 



v 



^ exp(— 7)r i( ^ 

2py(v + i) p 11 2 



-a 



+ 1 



, V+1 * 4p 2 > 



where I v (•) is the modified Bessel function of the first 
kind of order v, p’’ ( • ) is the Gamma function, and 
|F, (X, ,X 2 ;X S ) is the Confluent Hypergeometric function. In 
the special case of X ( =0, the . Confluent Hypergoemetric 

function ,F, (OjX^jX^) = 1 so that the integral in question 
can be simplified in this case. However, since the term S u 
includes the integration factor A, we may be able to calcu- 
late and express Equation 4 . 1 in a simpler form using the 
formula above for the case in which the variable A could be 
separated out in the term ^ Sn for a given jammer waveform. 

The same arguments apply to obtaining the probability of 
error under the hypotheses H 0 . Thus Pe can be shown to be 
given by 




fET x 2 

Q ( ^i£,x)x exp ( — y ) 



(4.2) 



Io (x ^-) dA 



dx 



dcp 
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Using Equations 4.1 and 4.2, we can evaluate the performance 
of receiver in a Rayleigh fading environment from 

= 4 (P el + P eo) 

assuming that the two signals are equally likely to be sent. 
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V. DESCRIPTION OF GRAPHICAL RESULTS 



A. GENERAL 

In this chapter, the analytical results of the previous 
chapters are now presented via graphical means based on the 
derived mathematical expressions for receiver probability of 
error . 

The plots presented display the receiver probability of 
error (Pe) as a function of SNR for the various jammer wave- 
forms previously considered for a set of JSR values. 

In each plot, the case JSR = 0 has been included in 
order to provide the basis for comparisons of the jammer 
effectiveness on the receiver performamce as it relates to 
additive white Gaussian noise only interference. 

B. ASK (ON - OFF KEYING) 

The graphical results for the incoherent ASK receiver 
performamce are presented first. These plots correspond to 
numerical evaluations of Equations 2.5 and 2.7 for equally 
likely hypothesis, that is, PCH,) and P(H 0 ) are equal to 
1 / 2 . 

The plot of Pe for ASK modulation is shown in Figure 6.7 
as a function of SNR for fixed values of JSR, using a jammer 
as specified in Equation 2.6 . Figure 6.7 clearly shows the 
'break point' phenomenon in which if JSR increases beyond a 
certain value ( 0.25 in this figure), Pe increases with 
increasing SNR. From this figure, one can observe that 16.0 
db of SNR is required to obtain a Pe of po -5 without jammer, 
i.e., at a JSR value of 0.0. In comparison, it takes 23.5 
db of SNR to obtain the same Pe for a JSR value of 0.1. 
Thus, in the presence of a jammer with a JSR value below the 
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break point, we need a larger SNR in order to obtain the 
same performance level of a receiver operating without a 
jammer interference. However, in a jamming environment, a 
JSR value above the break point produces a Pe which 
increases with increasing SNR. In fact, there is no value 
of SNR that can produce Pe of 10 -5 for a JSR above the break 
point. For the case of ASK modulation, the break point 
occurs at a value of JSR which is approximately 0.25, as 
obtained from Equations 2.8 and 2.9. 

Figure 6.8 shows the variable thresholding effect on the 
jamming situation with JSR = 0.3 beyond the break point (JSR 
= 0.25). Instead of using the fixed threshold as given by 
Equation 2.5, the variation of the receiver threshold 
obtained from Equation 2.5 by changing the value of R , can 
reduce the jamming effect over a restricted range of SNR 
values as shown. However, since as shown in Equation 2.13 
the variation of the value R does not significantly affect 
the threshold value, the variation over a wide range of 
values of R does not result in a significant change in Pe . 

C. FSK WITH TONE JAMMER 

This section presents graphical results pertaining to 
jamming effects on FSK modulation with a single tone jammer 
acting against one of the two channels and a jammer 
consisting of two different tones acting against both chan- 
nels simultaneously. 

Figure 6.9 corresponds to the performance of the optimum 
FSK receiver in which ^4 = 1/2. Equation 3.10 or Equation 
3.12 is used to evaluate performance with the near optimum 
jammer specified in Equation 3.11 . This jammer waveform 
can be thought of as 'mark' and 'space' channel jamming. 
Figure 6.9 shows a similar result to that found in the ASK 
case except that the breakpoint occurs at a higher value of 
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JSR than that found for ASK. This breakpoint occurs at a 
JSR somewhere between 0.5 and 1.0 as shown in this plot. 
From this figure it can be noted that 13.5 db of SNR is need 
to obtain a Pe of io -5 f° r a JSR value of 0.0, but the same 
Pe is obtained by increasing the SNR to 16.5 db for a JSR of 
0.1. This demonstrates that relatively low JSR values 
require a significant SNR boost in order to maintain a 
certain desired Pe value. 

As shown in Figures 6.7 and 6.9, comparison of ASK and 

FSK modulation reveals that FSK is somewhat less vulnerable 

to jamming. However it must be remembered that the jammer 

waveform n.(t) used in each case is different, 
o. 

The effect of the single tone jammer on the optimum FSK 
receiver is presented in Figure 6.10 which corresponds to 
the evaluation of Equation 3.15 . The single channel 
jamming on either the 'mark' or the 'space' channel has the 
same effect insofar as single tone jamming is concerned. 
Therefore the effect of 'mark' channel jamming only is eval- 
uated and plotted. Note that in Figure 6.10, 19.5 db of SNR 
is required to obtain a Pe of io -5 as compared to an SNR of 
24.5 db in Figure 6.9 for a JSR value of 0.3 with Pe = IQ -5 
also . 

As expected, Figures 9 and 10 demonstrate the fact that 
single channel jamming is less effective than simultaneous 
jamming of 'space' and 'mark' channels with a near optimum 
j ammer . 

The effect of a variable threshold on FSK will be 
considered by changing the value of q 4 in the modified 
receiver shown in Figure 6.6 . For a value of <yt, other than 
1/2, the simultaneous jamming of 'mark' and 'space' channel 
results in a compensation of the other channel such that the 
jamming effect remains the same as in the case of qC = 1/2. 
In other words, it is difficult to reduce the near optimum 
jamming effect by means of a varying the threshold. 



57 



On the other hand, when one channel jamming is applied 
to the receiver, the jamming effect can be reduced by 
adjusting the threshold with increasing SNR as shown in 
Figure 6.11, particularily for JSR of 1.0 and various values 
of A moment's reflection will reveal that the 'mark' 
channel jammer increases the output power level of the upper 
envelope detector (see Figure 6.6) so that when the 'space' 
signal is sent, the error increases. Therefore, for this 
type of error, under the assumption that the 'space' signal 
has been transmitted, Pe can be reduced by lowering the 
level of the output of the multiplier by using an appro- 
priate value of qC • 

D. FSK WITH FM JAMMER 

This section presents the effect of an FM sweep jammer 
using sinusodal modulation on noncoherent FSK signaling. The 
FM jammer was designed to sweep the bandwidth occupied by 
the signal several times during a bit interval. Thus the 
jammer effectiveness was investigated as a function of the 
number of times the jammer sweeps over the band of the 
signal during one bit-time interval. 

Figure 6.12 shows the result for one sweep of the jammer 
per bit interval. Figure 6.13 shows the result of 
increasing the sweeping to two sweeps per bit interval. 
These plots show that the FM sweep jammer can be more effec- 
tive by increasing the number of times of jammer sweeping 
during a bit interval. This can be expected from the results 
obtained in Equation 3.27 . In Figure 6.12, in order to 
obtain a Pe of 10 -5 , 16.0 db of SNR is required for JSR value 
of 0.3, but in Figure 6.13 which corresponds to twice the 
jammer sweeping, the same Pe can be obtained by increasing 
the SNR value to 19.5 db for the same JSR value of 0.3. In 
comparison to the previous case of FSK with tone jamming, 
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for the same jamming environment (i.e., JSR value of 0.3) 
Figures 6.9 and 6.10 show that 24.5 db and 19.5 db of SNR 
are required respectively in order to get the same Pe of 

10" 5 . 

These different requirements of SNR value for various 
jammer waveforms show that the FM sweep jammer can be effec- 
tive but in general is not as effective as the near optimum 
jammer. Note that from a practical point of view, the added 
complexity of FM jammer waveform may make it an unlikely 
candidate for a replacement of the near optimum jammer. 
However, one advantage the FM sweep jammer has over the near 
optimum jammer is that the former can spread its power over 
a large bandwidth easily and therefore is more effective 
than the latter in the case of lack of exact information or 
knowledge about the signal carrier frequency. 



59 



VI. CONCLUSIONS 



The familiar model in which Gaussian noise is the total 
interference is not adequate when jamming or interference 
signals are present in the transmission environment. This 
thesis has analyzed the effect of various deterministic 
jammer waveforms in terms of probability of error on binary 
incoherent receivers operating in the presence of noise. 

From the jammer point of view, the goal is to cause the 
maximum possible error to the various receivers while making 
efficient use of its available power (i.e., with fixed 
jammer power). For coherent receivers, it was proved that 
the optimum jammer waveform is made of a deterministic 
signal proportional to the difference of the binary signals 
used to carry the digital information. This thesis has 
demonstrated that those optimum jammers derived for coherent 
receivers perform their function as near optimum jammers 
satisfactorily against incoherent receivers. 

Therefore, these near optimum jammers can be concluded 
to be one of the most attractive candidates for efficient 
jamming of binary incoherent communication systems. An 
optimum jammer has not been derived or analyzed because the 
complexity of the expression for Pe makes it very difficult 
if not impossible to derive the optimum jammer waveform in 
closed form. Other jamming waveforms such as single channel 
jamming and FM sweep jamming showed its inferiority in 
comparison to near optimum jammer waveforms and their effi- 
ciency can be reduced by means of appropriate variations of 
the receiver threshold. 
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Figure 6.1 Correlation Receiver for Binary Signals. 
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Figure 6.2 Quadrature Receiver for ASK. 




Figure 6.3 Alternate Form of Quadrature Receiver. 
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Figure 6.4 Incoherent Matched Filter Receiver. 
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Figure 6.5' Incoherent Frequency Shift Keying (FSK) Receiver. 
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Figure 6.6 Modified Incoherent FSK Receiver. 
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Figure 6,7 Near Optimum Jammer for ASK Modulation, 
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Figure 6.8 Variable Threshold Effect on ASK Modulation. 
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Figure 6.9 Near Optimum Jammer for FSK Modulation. 
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Figure 6.10 Single Channel Jamming for FSK Modulation. 
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Figure 6.11 Variable Threshold Effect on Single Channel Jamming for FSK. 
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Figure 6.12 Frequency Modulated Jammer for FSK Modulation (K 
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Figure 6.13 Frequency Modulated Jammer for FSK Modulation (K 



APPENDIX A 

DIGITAL COMPUTER IMPLEMENTATION OF THE MARCUM Q- FUNCTION 



The Marcum Q- function occurs frequently in communication 
problems involving incoherent detection of signals with 
single or multiple observations in the presence of noise and 
jamming. So it is often necessary to compute values for the 
Marcum Q-function which is defined by 



Q( a , B 



2 2 

V exp(- - — ) I Q (a V) dV‘ 

2 

CP w 2 4. 2 

1 - V exp( - ) I 0 (a V) dV 

Jo 2 



(A. 1) 



where 1^ ( • ) is the modified Bessel function of zero order. 

It is noted that the integrand f ( V , pO of Equation A.l is 

♦ 

the Rician density function which is sometimes called the 
generalized Rayleigh density function. The normalized 
Rician distribution for y 2 - 1 is shown in Figure A.l . 

For large values of q( , the Rician density function f(V) 
can be approximated by the normalized Gaussian density func- 
tion with mean value of approximately ^ . This can be Justi- 
fied as follows 

2 2 

f(V) = V exp(- ) I Q ( a V) dV 
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Figvire A . 1 



Rician Density Function. 

_ x 



by use of the approximations I (X) 



[Z7 C>< ’ 



and = V 1 



for large ^ . 

In order to implement the computations of the Marcum 
Q- function using a digital computer, particular interest 
must be focused on values of the function on the tails of 
the Rician density defined by the integrand in Equation A.l 
. Specif ical ly , small values around the tails of the 
density must be monitored because of the limitations of the 
digital computer. By a change of variables (X = ^ V) the 
Marcum Q- function defined above can be expressed as 



Q( a , B ) = 1 - -K 

' a 

where f(X) -Q 1 Q (X) . Since the digital computer can treat 
the values in the range of y such that -180.218 5 y £ 

174.643 for the function e V ( for the IBM 3033 ), for 
digital computer implementation of the Marcum Q-function the 
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integral limits need to be adjusted to meet the conditions 
for these acceptable values in the computer. Negligible 
areas around the tails of the density beyond certain bounds 
must be discarded without significantly increasing overall 
error. Therefore from Equation A. 2, the integral limits can 
be substituted with appropriate bounds such that exponent 
satisfies the condition given by 



/ x , . 2 

- L eT . ~ ^ > -18 0 

In other words, the variable X must be located within the 

range pC 1 ~ ol /360 < X S ^ 1 * oC y 3 60. For a value of X 

beyond that range the exponential of Equation A. 2 is too 

-l?o 

small (less than Q ) compared to other computed values 
so that the area outside of the new limits can be neglected. 

Let us denote UL and LL the upper and lower limit 
respectively defined by 



UL 



2 

a 



LL = a 



2 



a/360 



Then we can consider the various situations case by case. 
If the lower limit LL is positive for certain values of £»( , 
then in the case of < LL , that is, £3 < (/. - /* 3 6 0 , the 

integral value can be ignored so that the value of ) 

can be assigned to be one. In the case of LL < S UL, 

that is, qC - J~2 60* < £3 < pC + y/ 3 6 0 , the Marcum Q- function 

can be approximated by 



Q( a , 6 ) = 1 




exp 




f(x) 



dx 



In the case of ^ > UL, that is, ^ > oi ♦ J 360, 
Marcum Q- function can be computed from 
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On the other hand if the lower limit LL is negative for some 
value of , then there are two cases to be cons idered. In 
the case of ^ ^ UL, that is, ^ < o( + / 360, the Marcum 
Q-function can be substituted with th^ value given by 
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ction can be computed 


r 


,x v 2 i 

(— - a) 


/ x exp 


a 


Jo 


L 2 J 



f(x) dx 

> </. ♦ ~/~ 360 ', 
from 

i 

f ( x ) dx 



Here it has been assumed that the function f(X) defined by 

m ft 

e I 0 (X) does not impose the limitation of computation on 
the digital computer and that the library functions for the 
computation of e Ij(X) and its integration with desired 
accuracy are available to the user. 
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